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A  series  of  1: 20-scale  model  tests  were  conducted  to  investigate  the  effects  of  tunnel 
diameter,  tunnel  length,  and  chamber  loading  density  on  the  Airblast  Inhabited 
Building  Distance  (IBD)  from  accidental  explosions  in  underground  magazines.  The 
results  indicate  that  the  IBD  decreases,  to  different  degrees,  with  increased  tunnel 
diameter,  decreased  loading  density,  or  increased  tunnel  length.  The  analysis  shows 
that  these  effects  are  primarily  a  function  of  the  total  volume  of  the  underground 
system. 


INTRODUCTION 

The  U.S.  Army  Engineer  Waterways  Experiment  Station  (WES)  has  conducted  research  on 
airblast  propagation  inside  underground  ammunition  storage  facilities.  The  purpose  of  the 
research  was  to  develop  magazine  design  enhancements  which  would  minimize  the  Airblast 
Inhabited  Building  Distance  (IBD)  should  an  accidental  internal  detonation  occur.  This  effort 
was  a  part  of  the  Joint  U.S. /Korea  R&D  Study  for  New  Underground  Ammunition  Storage 
Technologies. 

A  series  of  l:20-scale  experiments  were  conducted  at  WES  to  establish  a  baseline  data  base 
against  which  to  evaluate  the  effects  of  blast  reduction  designs  on  the  IBD.  Parameters 
evaluated  during  the  study  included  tunnel  length,  tunnel/chamber  diameter  ratio,  and 
chamber  loading  density  (Table  1).  The  parameter  study  was  performed  with  a  simple,  “shot 
gun  design  of  an  underground  magazine  (storage  chamber  and  access  tunnel  coaxially 
aligned). 

DESCRIPTION  OF  THE  MODEL 

The  small-scale  model  consisted  of  a  detonation  chamber  and  exit  tunnel.  The  detonation 
chamber  was  fabricated  from  a  50.8  cm  (inside  diameter)  steel  pipe  and  rings  of  steel  plate 
(7.6  cm  thick),  and  had  an  outer  diameter  of  81.3  cm.  The  exit  tunnel  was  made  of  heavy- 
walled  steel  pipe,  with  inside  diameters  0.3,  0.5,  and  0.7  times  the  chamber  diameter. 

Explosive  charges  for  the  model  tests  were  assembled  from  blocks  of  Composition  C-4. 
Nominal  dimensions  of  the  basic  charge  blocks  were  5.08  cm  wide,  by  2.54  cm  thick,  by  22.5 
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cm  long,  and  weighed  0.609  kg.  The  explosive  blocks  were  stacked  to  provide  explosive 
loading  densities  (charge  weight  per  cubic  meter  of  chamber  volume)  of  1.67,  3.33,  and  5.00 
kg/m3.  The  charges  were  positioned  at  the  center  of  the  chamber,  supported  on  a  light  wire 
frame.  An  RP-83  electric  bridge  wire  (EBW)  detonator  embedded  in  the  end  of  the  C-4 
charge  (end  closest  to  the  exit  tunnel)  initiated  the  explosive  charge. 

Instrumentation  to  measure  the  airblast  pressures  consisted  of  airblast  gages  located  along  the 
interior  surface  of  the  detonation  chamber  and  exit  tunnel,  and  external  airblast  pressure  gages 
located  on  the  ground  surface  outside  the  portal,  along  the  extended  centerline  of  the  exit 
tunnel.  Kulite  Corporation  Model  HKS-375  pressure  transducers  were  used  for  the  internal 
airblast  measurements.  The  internal  gages  were  flush-mounted  on  the  inner  wall  of  the 
detonation  chamber  and  pipe  sections  (side-on  pressure  measurements)  .  These  gages  had 
dynamic  ranges  of  210  MPa  to  3.4  MPa,  with  associated  resonant  frequencies  of  725  to  650 
kHz. 

The  external  measurements  employed  Kulite  Model  XT-190  airblast  gages.  These  instruments 
had  dynamic  ranges  between  13.8  MPa  and  35  kPa,  and  associated  resonant  frequencies 
between  650  kHz  and  70  kHz.  The  free-field  gages  were  mounted  flush  to  the  ground  surface 
beyond  the  tunnel  exits  at  ranges  (from  the  exit)  of  1.8,  3.0,  4.8,  and  8.0  m. 

BLASTX3  CALCULATIONS 

A  computational  model  of  the  small-scale  experiments  was  developed  using  BLASTX3,  a 
WES-developed,  empirical-based  airblast  code.  This  code  has  undergone  a  recent  revision  to 
allow  calculation  of  airblast  propagation  in  tunnels  for  moderate  to  large  loading  densities. 
Data  from  the  WES  small-scale  model  test  program  and  one-dimensional  hydrocode 
calculations  were  used  in  developing  the  improved  version  of  the  code  (BLASTX3). 

Two  experiments  (Series  1,  Tests  9  and  21),  both  with  tunnel/chamber  diameter  ratio's  of  0.5 
and  chamber  explosive  loading  densities  of  5  kg/m3,  were  calculated  with  the  BLASTX3 
computational  model.  The  length  of  the  exit  tunnels  were  2  m  (Test  9)  and  4  m  (Test  21). 

The  physical  model  was  divided  into  cylindrical  sections  with  an  approximate  length  to 
diameter  aspect  ratio  of  3:1  for  the  computational  model.  The  length-to-diameter  ratio  of  the 
model  detonation  chamber  was  3.5  to  1.  For  the  computation  model,  the  chamber  was 
modeled  as  a  single  cylindrical  room.  Similarly,  the  exit  tunnels  were  divided  into  3  (Test  9) 
and  6  (Test  21)  equal  cylindrical  sections  in  the  computational  model. 

Since  BLASTX3  can  only  model  spherical  charges,  the  rectangular  parallelepiped  explosive 
charge  used  in  the  physical  experiments  was  represented  by  a  single  spherical  charge  int  he 
computation. 


RESULTS  OF  EXPERIMENTS  AND  CALCULATIONS 


From  the  small-scale  experiments,  the  normalized  peak  overpressure  data  recorded  in  the 
free-field  beyond  the  tunnel  portal  are  plotted  versus  normalized  distance  in  Figures  2  through 
8.  The  non-dimensional  pressure  values  were  obtained  by  dividing  the  calculated  effective 
exit  pressure  by  the  measured  peak  free-field  pressure  (DDESB,  1993).  The  effective  exit 
pressure  (Pw)  is  calculated  from 


EQUATION 


Pw  =  1770.5 


1,35 

3 


where 

Q  is  the  net  explosive  weight  in  the  chamber,  kg,  and  Vt  is  the  total  volume  of  the 
underground  storage  system,  m3. 

The  normalized  distances  were  computed  by  dividing  the  distances  from  the  exit  portal  by  the 
hydraulic  tunnel  diameter  for  turbulent  flow,  which  is  obtained  from 


D  = 


4A 

T 


where 

A  is  the  minimum  cross-sectional  area  of  the  tunnel,  m2,  and  P  is  the  perimeter  of  the 
minimum  tunnel  cross-section,  m. 

The  equation  developed  by  the  DOD  Explosives  Safety  Board  (DDESB)  to  predict  the 
normalized  IBD  is  included  in  Figures  2  through  8  to  provide  a  reference  against  which  to 
evaluate  the  experimental  data.  The  normalized  IBD  prediction  relation  is 

^  =  (R/D)135 


where 


P  is  the  peak  free-field  overpressure,  kPa,  and  R  is  the  distance  from  the  portal,  m. 


Although  there  is  some  data  scatter,  the  normalized  free-field  pressure  shows  no  significant 
influence  of  the  chamber  loading  density  (Figures  2,  3,  and  4,  respectively.  As  shown  in 
these  graphs,  the  trend  of  the  experimental  data  follows  the  IBD  prediction  line.  The 
normalized  peak  pressures  in  the  free-field  as  calculated  by  BbASTX3  are  also  shown  in 
Figure  3,  (for  the  5.00  kg/m3  loading  density).  The  BLASTX3  curve  follows  the  trend  of  the 
data  (triangles)  for  this  loading  density. 

The  effect  of  tunnel  length  on  normalized  free-field  overpressure  is  shown  in  Figures  5,  6,  and 
7.  Although  the  experimental  results  exhibit  some  scatter,  no  significant  influence  of  tunnel 
length  can  be  identified  from  these  test  results.  Figure  6  shows  the  normalized  overpressures 
calculated  by  BThASTX3,  for  2  an  4-m  tunnel  lengths.  The  BLASTX3  results  indicate 
decreased  free-field  pressures  with  increased  tunnel  lengths. 

A  comparison  of  data  from  the  three  tunnel/chamber  diameter  ratios  that  were  evaluated  in  the 
small-scale  test  program  is  shown  in  Figure  8.  The  data  exhibits  some  scatter,  masking  any 
effects  of  tunnel  diameter.  Base  on  these  results,  the  effect  of  tunnel  diameter  on  free-field 
pressures  (and  the  IBD)  was  judged  to  be  insignificant. 

As  mentioned  earlier,  the  total  volume  of  the  underground  magazine  system  was  used  to 
calculate  the  effective  exit  pressure.  In  effect,  normalizing  the  free-f  ield  pressure  data  with 
the  effective  exit  pressure  accounts  for  the  extra  volumes  associated  with  larger  or  longer 
access  tunnels.  Similarly,  the  effective  exit  pressure  calculation  includes  the  charge  weight  or 
mass.  Thus,  the  normalizing  procedure  essentially  eliminates  the  effect  of  loading  density 
from  the  data  comparisons  presented  in  this  paper.  Therefore,  a  doubling  of  the  explosive 
load  in  an  underground  magazine  should  increase  the  free-field  pressures  by  a  factor  of  two 
taken  to  the  0.45  power. 

CONCLUSIONS 

The  results  of  explosive  tests  conducted  in  l:20-scale  models  provided  indications  of  the 
extent  to  which  the  IBD  decreases  as  a  function  of  increased  tunnel  diameter  and  tunnel 
length,  and  decreased  loading  density.  However,  when  the  peak  airblast  pressures  are 
normalized  and  plotted  against  normalized  distance  from  the  tunnel  portal  (using  current 
DDESB  formulas  for  normalizing  these  parameters),  there  is  no  significant,  independent 
effect  of  chamber  loading  density,  tunnel  length,  or  tunnel  diameter.  These  conclusions  were 
supported  by  the  results  of  calculations  using  the  BLASTX3  computer  model. 
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TABLE  1.  JOINT  U.S./ROK  R&D  PROGRAM  FOR  NEW  UNDERGROUND 
AMMUNITION  STORAGE  TECHNOLOGIES:  SMALL-SCALE 
MODEL  TEST  PARAMETERS 


TEST 

NO. 


VENT  PIPE 
LENGTH 
(m) 


TUNNEL / CHAMBER 
DIAMETER 
RATIO 


LOADING  DENSITY 


(kg/m3) 


TABLE  1.  JOINT  U.S./ROK  R&D  PROGRAM  FOR  NEW 
UNDERGROUND  AMMUNITION  STORAGE  TECHNOLOGIES: 
SMALL-SCALE  MODEL  TEST  PARAMETERS 


FIGURE  1. 

CONTROL  MAGAZINE  TEST  SERIES,  PLAN  VIEW  OF  GENERIC 
TEST  LAYOUT.  VENT  PIPE  LENGTHS  RANGED  FROM  1  TO  4  M. 
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Note:  *  denotes  stagnation  pressure  measurement 
gage  placed  surface  flush  1 .3  cm  in  front 
of  2.5  cm  cube 


Figure  1 .  Control  magazine  test  series,  plan  view  of  generic  test  layout. 
Vent  pipe  lengths  ranged  from  1  to  4  m. 
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Figure  2.  Effect  of  loading  density:  comparison  of  normalized  peak  free- 
field  overpressure  from  small-scale  model  tests  at  chamber  loading 
densities  of  1.67,  3.33,  and  5.00  kg/m3  for  a  tunnel/chamber  diameter  ratio 

of  0.3,  and  tunnel  length  of  4  m. 


Ld 

C£ 

Z) 

00 

in 

LlI 

Q_ 

< 

Ld 

Q_ 

Q 

Ld 

ISI 


< 


o 


Figure  2.  Effect  of  loading  density:  comparison  of  normalized 
peak  free-field  overpressure  from  small-scale  model 
tests  at  chamber  loading  densities  of  1.67,  3.33,  and 
5.00  kg/m3  for  a  tunnel/chamber  diameter  ratio  of  0.3, 
and  tunnel  length  of  4  m. 


Figure  3.  Effect  of  loading  density:  comparison  of  normalized  peak  free- 
field  overpressure  from  small-scale  model  tests  at  chamber  loading 
densities  of  1.67,  3.33,  and  5.00  kg/m3  for  a  tunnel/chamber  diameter  ratio 

of  0.5,  and  tunnel  length  of  4  m. 
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Figure  3.  Effect  of  loading  density:  comparison  of  normalized 
peak  free-field  overpressure  from  small-scale  model 
tests  at  chamber  loading  densities  of  1.67,  3.33,  and 
5.00  kg/m3  for  a  tunnel/chamber  diameter  ratio  of  0.5, 
and  tunnel  length  of  4  m. 


NORMALIZED  PEAK  PRESSURE 


Figure  4.  Effect  of  loading  density:  comparison  of  normalized  peak  free- 
field  overpressure  from  small-scale  model  tests  at  chamber  loading 
densities  of  1.67,  3.33,  and  5.00  kg/m3  for  a  tunnel/chamber  diameter  ratio 

of  0.7,  and  tunnel  length  of  4  m. 


Figure  4.  Effect  of  loading  density:  comparison  of  normalized 
peak  free-field  overpressure  from  small-scale  model 
tests  at  chamber  loading  densities  of  1.67,  3.33,  and 
5.00  kg/m3  for  a  tunnel/chamber  diameter  ratio  of  0.7, 
and  tunnel  length  of  4  m. 


NORMALIZED  PEAK  PRESSURE 


Figure  5.  Effect  of  tunnel  length:  comparison  of  normalized  peak  free- 
field  overpressure  from  small-scale  model  tests  with  access  tunnel  lengths 
of  1,  2,  and  4  m  for  a  chamber  loading  density  of  5.00  kg/m3  and  a 
tunnel/chamber  diameter  ratio  of  0.3. 


Figure  5.  Effect  of  tunnel  length:  comparison  of  normalized  peak 
free-field  overpressure  from  small-scale  model  tests 
with  access  tunnel  lengths  of  1,  2,  and  4  m  for  a 
chamber  loading  density  of  5.00  kg/m3  and  a 
tunnel/chamber  diameter  ratio  of  0.3. 


Figure  6.  Effect  of  tunnel  length:  comparison  of  normalized  peak  free-field 
overpressure  from  small-scale  model  tests  with  access  tunnel  lengths  of  1, 
2,  and  4  m  for  a  chamber  loading  density  of  5.00  kg/m3  and  a 
tunnel/chamber  diameter  ratio  of  0.5. 
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NORMALIZED  DISTANCE  FROM  PORTAL 

Figure  6.  Effect,  of  tunnel  length:  comparison  of  normalized  peak 
free-field  overpressure  from  small-scale  model  tests 
with  access  tunnel  lengths  of  1,  2,  and  4  m  for  a 
chamber  loading  density  of  5.00  kg/m3  and  a 
tunnel/chamber  diameter  ratio  of  0.5. 


NORMALIZED  PEAK  PRESSURE 


Figure  7.  Effect  of  tunnel  length:  comparison  of  normalized  peak  free-f 
ield  overpressure  from  small-scale  model  tests  with  access  tunnel  lengths  of 
11  2,  and  4  m  for  a  chamber  loading  density  of  5.00  kg/m3  and  a 
tunnel/chamber  diameter  ratio  of  0.7. 


Figure  7.  Effect  of  tunnel  length:  comparison  of  normalized  peak 
free-field  overpressure  from  small-scale  model  tests 
with  access  tunnel  lengths  of  1,  2,  and  4  m  for  a 
chamber  loading  density  of  5.00  kg/m3  and  a 
tunnel /chamber  diameter  ratio  of  0.7. 


Figure  8.  Effect  of  tunnel/chamber  diameter  ratio:  comparison  of 
normalized  peak  free-field  overpressure  from  smallscale  tests  with 
tunnel/chamber  diameter  ratio's  of  0.3,  0.5,  and  0.7  for  a  chamber  loading 
density  of  5.00  kg/m3  and  tunnel  length  of  4  m. 
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Figure  8 


Overpressure  generated  (in  psi)  when  sandbags  are  used. 


